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GOALS  &  OBJECTIVES 

Turbulence  in  the  surfzone  and  nearshore  mixes  momentum  vertically,  transmits  stress  to  the  sea 
bed,  influences  the  structure  of  the  cross-  and  alongshore  currents,  and  controls  the  suspension  of 
sediment  from  the  sea  bed.  In  many  coastal  and  shelf  environments,  the  sea-bed  is  the  primary 
source  of  turbulence  due  to  bottom  induced  shear.  In  the  surfzone,  the  breaking-wave  generated 
turbulence  likely  dominates  over  bottom  generated  turbulence.  However,  the  dynamics  of  turbu¬ 
lence  in  the  nearshore  and  surfzone  under  breaking  waves  is  poorly  understood. 

Realistic  three-dimensional  simulations  of  surfzone  hydrodynamics  and  sediment  transport, 
which  are  currently  being  attempted,  for  example,  in  the  recent  nearshore  NOPP  project,  will  not 
be  possible  without  at  least  a  rudimentary  understanding  of  nearshore  turbulence  dynamics.  Long 
term  goals  include  addressing  some  of  the  unresolved  science  issues  through  analysis  and  modeling 
of  existing  field  measurements  to  quantify  turbulence  dynamics. 

To  address  nearshore/surfzone  turbulence  dynamics,  this  research  had  two  approaches,  model¬ 
ing  work  and  analysis  of  field  data,  which  has  resulted  in  3  journal  publications  {Feddersen  and 
Trowbridge,  2005;  Feddersen  and  Williams,  2006;  Feddersen  et  ah,  2006). 


APPROACH 

The  model  development  (programming)  and  testing  work  was  done  by  the  PI  (Feddersen),  The 
basis  of  the  coupled  turbulence-bubble  model  is  a  standard  k-t  model  augmented  with  a  wave¬ 
breaking  turbulence  source  terms  in  an  approach  similar  to  that  used  for  open-ocean  wave  breaking 
(Craig  and  Banner,  1994),  but  with  the  time-dependence  of  breaking  retained.  The  model  is 
described  in  detail  in  Feddersen  and  Trowbridge  (2005). 

The  field  campaign  consisted  of  a  main  instrumented  frame  was  deployed  for  2  weeks  at  the 
U.S.  Army  Corps  of  Engineers  Field  Research  Facility  (FRF)  at  Duck  N.C.  in  approximately  3.2  m 
mean  water  depth  with  the  variable  tides  and  wave  heights  (Fig.  1).  At  this  depth,  there  was  white¬ 
capping  wave  breaking  but  no  depth-limited  wave  breaking.  The  main  frame  was  instrumented 
with  a  vertical  array  spanning  most  of  the  water  column  of  3  ADV  (Acoustic  Doppler  Velocimeter) 
to  to  measure  velocity.  These  measurements  of  the  vertical  structure  of  turbulence  in  the  nearshore 
are  unique  and  had  not  before  been  attempted  in  the  nearshore  region. 

Field  data  analysis  included  estimating  the  dissipation  from  the  ADVs,  using  the  frozen  turbu¬ 
lence  hypothesis  and  the  inertial-dissipation  technique  (Trowbridge  and  Elgar,  2001).  The  verti¬ 
cally  separated  ADVs  were  used  to  infer  the  turbulent  momentum  fluxes,  Le.,  v'w'  (e.g.,  Feddersen 
and  Williams,  2006)  that  are  part  of  the  shear  production  term  in  the  TKE  dynamics. 
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Figure  1 :  Schematic  of  the  ADV  locations.  The  view  is  toward  offshore  (+a:),  and  the  vertical  ^ 
and  alongshore  y  coordinates  are  indicated.  ADV  1  is  upward  looking.  The  vertical  location  of  the 
ADV  sensing  volumes  (indicated  by  the  small  circle)  are  given 


RESULTS 

This  grant  has  resulted  in  3  journal  publications,  detailed  individually  below. 


Surfzone  TUrbulence  Modeling 

The  newly  surfzone/nearshore  turbulence  model  is  able  to  reproduce  the  observed  non-dimensionalized 
turbulence  dissipation  c/{g^hy^‘^  from  three  different  data  sets.  The  nondimensionalization  of  dis¬ 
sipation  collapses  the  three  data  sets  which  come  from  different  beaches  and  wave  conditions 
{Feddersen  and  Trowbridge,  2005). 


Estimating  Reynolds  Stresses 

An  important  component  of  testing  a  turbulence  model  is  developing  field  estimates  of  the  Reynolds 
stress  component  v'w'.  As  is  discussed  in  Trowbridge  (\99S),  estimation  of  v'w'  in  surface  gravity 
wave  dominated  environments  requires  special  methods.  Trowbridge  {\99&)  developed  a  v'w'  esti¬ 
mation  method  using  two  current  meters  and  tested  it  near  the  bed  with  two  horizontally  separated 
current  meters  {Trowbridge,  1998;  Trowbridge  and  Elgar,  2001).  Shaw  and  Trowbridge  (2001) 
developed  another  method  for  vertically  separated  sensors  with  weak  wave  conditions  near  the 
bed  off  of  the  continental  shelf.  However,  for  the  vertical  ADV  array  (Fig.  1),  a  new  v'w'  estima¬ 
tion  method  had  to  be  developed,  because  of  problems  with  the  other  methods.  The  Trowbridge 
(1998)  method  has  significant  wave  bias  problems  (not  shown).  The  Shaw  and  Trowbridge  (2001) 
method  v'w'  were  better  but  often  failed  a  test  of  the  integrated  velocity  cospectrum.  The  new 
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Figure  2:  Inter-sensor  comparison  of  Feddersen  and  Williams  (2006)  estimated  v'w'  for  all  combi¬ 
nations  of  vertical  location.  The  correlations  and  number  of  good  colocated  estimates  are  given. 


method,  described  in  Feddersen  and  Williams  (2006),  is  able  to  reproduce  a  constant  stress  layer 
(i.e.,  vertically  uniform  v'w')  in  the  wind-driven  nearshore  environment  (Fig.  2). 

This  work  Feddersen  and  Williams  (2006)  is  in  press  to  J.  Atmos.  Oceanic  Tech. 


Vertical  Structure  of  I\irbulent  Dissipation  in  the  Nearshore 

The  next  componenet  of  this  research  concerns  the  vertical  structure  of  turbulent  dissipation  in  the 
nearshore  and  is  now  in  press  to  JPO  {Feddersen  et  ai,  2006). 

The  dissipation  of  turbulent  kinetic  energy  (e)  is  estimated  with  a  new  method  {Feddersen 
et  al.,  2006)  at  each  of  the  ADVs  (Fig.  1).  The  mean  plus  the  first  EOF  captures  the  majority 
(91%)  of  the  e  variability.  The  observed  e  has  a  maximum  at  the  uppermost  ADV,  a  minimum 
at  the  middle  ADV,  and  another  maximum  but  less  than  ADV  3,  at  the  bottom  (Fig.  3).  This  is 
consistent  with  surface  generated  turbulence  due  to  wind-induced  wave-whitecapping  influencing 
ADV  3  and  near-bed  generated  turbulence  influencing  ADV  1  (see  for  example  Feddersen  and 
Trowbridge,  2005). 

Further  demonstrating  that  the  increased  surface  e  is  not  due  to  surfzone  depth-limited  wave 
breaking,  the  c  observations  are  scaled  with  a  surfzone  dissipation  scaling  {Feddersen  and  Trow¬ 
bridge,  2005)  giving  the  non-dimensional  c/{g^h)^^‘^.  At  the  same  relative  locations  in  the  water 
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Figure  3:  Vertical  structure  of  the  log-eof  e:  mean  (dashed)  and  mean  ±  first  EOF  (solid). 

column,  these  values  were  much  less  that  other  scaled  surfzone  c/{g^hy^^  {George  et  ai,  1994; 
Bryan  et  ai,  2003).  Another  process  is  giving  the  increased  near  surface  e. 

Wind-induced  wave  whitecapping  is  likely  that  process.  In  a  deep  lake,  Terray  et  al.  (1996) 
found  increased  e  near  the  surface  that  scaled  as 

e/4ig/(au^)  =  0.3(////,ig)-2  (1) 

where  z'  is  the  distance  from  the  surface,  u»  is  the  friction  velocity  and  a  is  a  constant.  This  im¬ 
plies  that  white-capping  was  inputing  turbulence  at  the  surface  which  diffuses  downward  resulting 
in  enhanced  e  relative  to  law  of  the  wall  scaling.  The  top  two  ADVs  follow  this  deep-water  scal¬ 
ing  (Fig.  4a,b).  This  demonstrates  that  the  processes  just  offshore  of  the  surfzone  that  result  in 
increased  near  surface  dissipation  are  the  same  wind-induced  white-capping  wave  breaking  as  in 
the  deep  water  {e.g.,  Terray  et  ai,  1996).  However  the  near  bed  ADV  does  not  follow  this  scaling 
and  is  likely  influenced  by  bottom  boundary  layer  processes. 
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